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The red blazar PMN J2345-1555 becomes blue 
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ABSTRACT 

The Flat Spectrum Radio Quasar PMN J2345-1555 is a bright 7-ray source, that recently 
underwent a flaring episode in the IR, UV and 7-ray bands. The flux changed quasi simul- 
taneously at different frequencies, suggesting that it was produced by a single population of 
emitting particles, hence by a single and well localized region of the jet. While the overall 
Spectral Energy Distribution (SED) before the flare was typical of powerful blazars (namely 
two broad humps peaking in the far IR and below 100 MeV bands, respectively), during the 
flare the peaks moved to the optical-UV and to energies larger than 1 GeV, to resemble low 
power BL Lac objects, even if the observed bolometric luminosity increased by more than 
one order of magnitude. We interpret this behavior as due to a change of the location of the 
emission region in the jet, from within the broad line region, to just outside. The correspond- 
ing decrease of the radiation energy density as seen in the comoving frame of the jet allowed 
the relativistic electrons to be accelerated to higher energies, and thus produce a "bluer" SED. 

Key words: BL Lacertae objects: general — quasars: general — radiation mechanisms: non- 
thermal — gamma-rays: theory — X-rays: general 



1 INTRODUCTION 

PMN 12345-1555 is a Flat Radio Specti-um Quasar (FSRQ) at 
z — 0.621 (Healey et al. 2008), bright in 7-rays, and detected by 
the Large Area Telescope (LAT) onboard the Fermi satellite in the 
first 3 months of its operation (Abdo et al. 2009), and later present 
in the ILAC (catalog of AGN detected in the first 11 months of 
operations, Abdo et al. 2010) and 2LAC (first 2 years, Nolan et al. 
2012). The 7-ray spectra were described by simple power laws. 
The fluxes above 100 MeV are reported in Tab.[T] together with the 
photon spectral indices. It can be seen that the source varied on long 
timescales, even if the reported fluxes are averaged over the entire 
period of observations. On the other hand the spectral index F was 
always greater than 2, indicating a i^Fi, peak at energies smaller 
than 100 MeV. This is typical for the FSRQ class (see Abdo et al. 
2010; Ghisellini, Maraschi & Tavecchio 2009). 

At the beginning of 2013 the source underwent a major flare, 
reaching, on Jan 13, 2013, a flux above 100 MeV of (110 ± 20) x 
10~* ph cm~^ s~^, more than 15 times the average photon flux 
listed in the 2LAC catalog (Tanaka 2013). Most interestingly, as 
noted in the telegram by Tanaka (2013), the spectral index became 
harder, F = 1.78 ± 0.13. This implies a rising spectrum in uF,^, 
which is instead typical of low power blazars, namely high energy 
peaked BL Lac objects. Since the spectrum hardened, the lumi- 
nosity variation between Jan 13 2013 and the average value in the 
2LAC catalog was almost a factor 30. The spectral change was ac- 
companied by a flaring behavior not only in 7-rays, but also in IR, 



Email: gabriele.ghisellini@brera.inaf.it 



Catalog/Date 


F_g [>0.1 GeV] 


Flat 


log Lj 


LBAS 


10.3 ± 1.3 


2.42 ±0.12 


46.89 


ILAC 


3.5 


2.37 ±0.10 


46.44 


2LAC 


6.5 


2.19 ±0.04 


46.79 


Jan 13 2013 


110 ± 20 


1.78 ±0.13 


48.25 



Table 1. Fluxes F_8 in units of 10"* ph cm-^ g-i above 100 MeV 
and coiTesponding spectral indices as reported in different catalogs, and in 
Tanaka et al. (2013) for the Jan 13, 2013 data. The last column repoits the 
K-coiTected [0.1-100 GeV] luminosities. 



optical and UV (Carrasco et al. 2013; Donate et al. 2013), where 
the source increased its flux by almost 2 orders of magnitude, and 
in X-rays (Donato et al. 2013), albeit with a more modest increase. 

Although other FSRQs with a relatively hard 7-ray spec- 
trum have been already detected by Fermi (discussed in Padovani, 
Giommi & Rau 2012 and Ghisellini et al. 2012), it is the first 
time that we can well document, in a given FSRQ (i.e. a blazar 
with relatively strong broad emission lines), the change from a soft 
(Flat > 2) to a hard (Flat < 2) 7-ray slope, accompanied by 
an intense flare of the source from the near IR to the 7-ray band. 
Previous hardening of the 7-ray spectrum have been observed in 
the FSRQ 4C +21.35 (=1222+216, Tanaka et al. 2011, interpreted 
by Foschini et al. 201 1 as a change in the location of the dissipation 
region in the same source) and in BL Lac (Bloom et al. 1997). 

The aim of this letter is to analyze the available data concern- 
ing the flare, to compare them to the existing older observations, 
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and to construct a coherent picture of what caused the "red to blue" 
transition of this blazar. 

We use a cosmology with h = SIa = 0.7 and SIm = 0.3. 



2 DATA ANALYSIS 

2.1 Swift observations 

Swift XRT and UVOT data were analyzed by using the HEASof t 
V. 5.13 software package with the CALDB updated on 21 Jan- 
uary 2013. XRT data were processed with xrtpipeline v. 
0.12.6 with standard parameters. The extracted count spectra 
have been grouped to have at least 25 counts per bin, in order 
to adopt the test and analyzed with xspec v. 12 . 8 . in 
the 0.3-10 keV energy band. The basic adopted model was a 
redshifted power law with a fixed Galactic absorption column of 
Nn = 1.64 X 10^*^ cm"^ (Kalberla et al. 2005). However, in the 
2013 Jan 14 observation (obsID 00038401009) a broken power law 
model is required at 99.53% according to the ftest. See Tab.|2] 

OpticalAJV photometry was done with UVOT, by adopting a 
5"-sized region for the source and a 7"~60" source-free annulus 
for the background. UVOT tasks uvotimsum and uvotsource 
were used to perform the analysis. The observed magnitudes, re- 
ported in Tab. |2] were then corrected for the Galactic absorption 
according to the filter-specific formulas of Cardelli et al. (1989) 
and converted into physical units flux by using the calibration of 
Poole et al. (2008) and Breeveld et al. (2010). 

2.2 Fermi/LAT observations 

Publicly available Fermi LAT data were retrieved from the Fermi 
Science Support Center (FSCC) and analyzed by means of the LAT 
Science Tools v. 9.27.1, together with the Instrument Response 
Function (IRF) Pass 7 and the corresponding isotropic and Galac- 
tic diffuse background models. Source (class 2) photons in the 0. 1- 
100 GeV energy range, collected on January 11-14 and coming 
from direction within 10° from the nominal source position were 
selected and filtered through standard quality cuts. Standard anal- 
ysis steps compliant with the FSSC recommendations were then 
performed. Besides the target and backgrounds, all the 2FGL point 
sources in the field were iwncluded in the model. PMNJ2345-1555 
was detected with high statistical significance (TS = 340, Mattox et 
al. 1996) and a high flux [Fo.i-iooGoV = (8.76 ± 0.16) x 10"^ 
photons cm~^ s"^] and a flat photon index F = —2.02 ± 0.01. A 
significant spectral break at _Eb = 2 ± 0.5 GeV is found fitting the 
data with a broken power law model, with low and high spectral in- 
dices Fi = 1.84 ± 0.15 and r2 = 2.58 ± 0.34. A similar analysis, 
but made in logarithmically spaced energy bins, gives results con- 
sistent with this break (solid black points in Fig.[T). This suggests 
that the high energy hump in the SED peaks at ~2 GeV. 



3 SPECTRAL ENERGY DISTRIBUTION AND MODEL 

Fig. [T] shows the SED of the source, and highlights the dramatic 
flux variations occurred from the IR to the 7-ray bands. The cyan 
bow-tie corresponds to the maximum seen by Fermi on Jan 13, 
2013; black symbols correspond to quasi simultaneous data taken 
between 1 1 and 14 Jan {Fermi), Jan 14 (X-rays and optical-UV by 
Swift). The brown symbols refer to Jan 17, 2013. The red symbols 
refer to two Swift observations merged together taken on Dec. 23, 
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17.66±0.05 


13.32±0.04 


13.71±0.04 
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17.78±0.04 


13.49±0.04 


13.88±0.04 



Table 2. Top: Results of the X-ray analysis. For all spectra a fixed Galac- 
tic A^H = 1-64 X 10^'^ cm~^ was assumed. The flux F0.2-10 is the 
unabsorbed flux in the [0.2-10 keV] band in units of lO"^-^ erg cm"-^ 
s^^. Bottom: UVOT Observed magnitudes, a: Average of 2 observations 
on 23/12/2008 and 10/01/2009. 
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Figure 1. The SED of PMN J2345-1555 and the fitting models. Black 
points correspond to Swift/UYOT and XRT data taken on Jan. 14, 2013 
and to Fermi data integrating the flux between Jan 1 1 and 14. The brown 
points are the UVOT data taken on Jan 17, 2013. The cyan bow-tie corre- 
sponds to the maximum seen by Fermi on Jan 13, 2013. We can compare the 
high state of Jan 2013 with the low state observed in Dec 23, 2008 by Swift 
(red points), while the red bow-tie refer to the average state of the first 3 
months of Fermi) operations (LBAS). Green symbols refer to archival data. 
We show the models for the two states, as explained in the text and with 
parameters listed in Tab. |3] with the different components labeled. 

2008 and Jan 10, 2009. The red bow tie in 7-rays is the average 
spectrum observed in the first 3 months of operations of Fermi. 

One of the most interesting feature of this SED concerns the 
X-ray spectrum: being a broken power law, it is strongly suggest- 
ing that in this band we see the end of the steep synchrotron compo- 
nent and the emergence of (presumably) the flat inverse Compton 
(IC) contribution. This is a defining property (Padovani &. Giommi 
1995) of intermediate BL Lacs (see Tagliaferri et al. 2000 for a 
similar spectrum in the BL Lac source ON 231). Since in the previ- 
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ous 5w(/j/XRT observation the entire X-ray spectrum was flat, and 
could be interpreted as IC, we must conclude that the peak of the 
synchrotron flux shifted towards higher frequencies. As a result, 
also the IC bump is "bluer", making the FennifLAT spectrum flat, 
and peaking at hi^c ~2 GeV. At the same time, the UVOT data are 
consistent with a synchrotron peak at i/s ~ 10^'' Hz. 

If the IC process with seed photons of frequency i^ext pro- 
duced externally dominates (External Compton, EC for short), we 
have (see Tavecchio & Ghisellini 2008): 
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where F is the bulk Lorentz factor; S = l/[r(l — /3cos6v]; 9^ 
is the viewing angle and 7pcakWicC^ is the energy of the electrons 
emitting at the peaks of the SED. Setting uc = 4.8 x 10^^ Hz (i.e. 
2 GeV) and T = 5, we derive 7poak ~ 700(r/15) or 7000(r/15) 
according if i^ext is the Hydrogen Lya or the peak of the IR emis- 
sion produced by the torus, assumed to be at 3 x 10^"^ Hz. The 
synchrotron peak frequency i^s ~ 10^^ Hz then yields: 



i^s = 3.6 X lo^^E^^ 
l + z 



B5 



Hz 



B > 



G (2) 



1.1 X lOVc 

This gives B ~ 70(r/15) or 0.7(r/15) G if the seed photons 
are coming from the broad line region (BLR) or from the IR torus, 
respectively. The very same relations hold for the SED observed 
when the 7-ray spectrum is steep, but with different i^s and vq . As 
we shall see, a coherent scenario explains the large shift in peak 
frequencies as due to the dissipation region, usually located within 
the BLR, moving beyond it. In such a case the reduced cooling 
allowed the electrons to reach larger 7pcak, corresponding to larger 
Us and vc- We then apply a model along these lines. 



3.1 The model 

We use the model described in detail in Ghisellini & Tavecchio 
(2009), and used in our previous blazar studies. For self consis- 
tency, and to explain the main parameters listed in Tab.|3] we repeat 
here the main properties of the model. 

The model accounts for several contributions to the radiation 
energy density, and how these and the magnetic one scale with the 
distance i?diss from the black hole of mass M. We consider radia- 
tion from the disk (i.e. Dermer & Schlickeiser 1993), the BLR (e.g. 
Sikora, Begelman & Rees 1994), a dusty torus (see Blazejowski et 
al. 2000; Sikora et al. 2002), the host galaxy light and the cosmic 
background radiation. It is a one-zone, leptonic model. 

The emitting region is spherical, of size R ~ O.l-Rdiss- The 

1/2 

jet accelerates in its inner parts with F oc ii^iss "^P '■° ^ value Fmax. 

The electron injection function Q{'y) [cm"'' s~^] is assumed 
to be a smoothly joining broken power-law, with a slope (3(7) oc 
7^°^ and 7^°^ below and above a break energy 7b: The total 
power injected into the source in relativistic electrons is P/ = 
mecV / Q(7)7d7, where V = (47r/3)i2^ 

The BLR, reprocessing 10% of the disk luminosity L^, is as- 
sumed to be a thin spherical shell located at a distance i?BLR = 
d^i The radiation energy density of the broad lines is 
constant within the BLR, but it is seen amplified by ~ F^ by 
the moving blob (as long as 7?diss < iiBLft). For illustration. 
Fig. [3] shows also the case of a "ring-like" BLR, lying in the 
disk, shaped as a torus of radius i?BLR and cross sectional radius 
equal to CI/Zblr- The (/blr profile is somewhat smoother in this 
case. The two cases (spherical and ring-like) should bracket the 
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Figure 2. Zoom of the SED of PMN J2345-1555 to better show the IR- 
UV part of the spectram. The 5w(/f/UVOT data of Dec. 23, 2008 can be 
described by a standard accretion disk model. Note the upturn in the optical 
spectrum taken by Shaw et al. (2009), labeled S09. To illustrate the possible 
uncertainties on the black hole mass, we show three accretion disk (+torus) 
models that have the same disk luminosity but different black hole masses, 
as labelled. In Jan 13, 2013 the optical-UV synchi'otron flux increased by 
almost 2 orders of magnitude, lai'gely overtaking the thermal emission. 



possible geometrical possibilities. A dusty torus, located at a dis- 
tance _RiR — 2.5 X 1Q^^L]/1^ cm, reprocesses ~30% of Ld in 
the far IR. In the inner parts of the accretion disk there is an X- 
ray emitting corona of luminosity Lx ~ 0.3Ld with a spectrum 
F{u) oc i/-^exp(-/ii//150keV). 

The energy densities of all these external components are cal- 
culated in the jet comoving frame, and used to calculate the result- 
ing EC spectrum. The internally produced synchrotron emission is 
used to calculate the synchrotron self Compton (SSC) flux. 

We adopt a standard Shakura & Sunjaev (1973) accretion disk 
spectrum. It depends on the black hole mass Al and the accretion 
rate M, that can be found if the SED shows signs of disk emission. 
In this case the total disk luminosity id fixes A/, and the peak fre- 
quency of the disk spectrum fixes M. Furthermore, we are helped 
by the presence of broad lines, that can be used as a proxy for Ld. 
This method, detailed in Calderone et al. (2013), returns very accu- 
rate black hole masses if the peak of the disk emission is visible. 



4 RESULTS 

Fig. [T] shows the result of the modeling: the red solid curve is the 
model for the data of Jan 13, 2013, and the solid blue line is the 
model fitting the Dec. 2008 UVOT and XRT data, together with 
the average 7-ray flux of the first 3 months of Fermi. Both models 
share the same black hole mass, accretion rate and viewing angle. 
All used parameters are listed in Tab.|3]that also reports, for ease of 
the reader, the set of parameters used in Ghisellini et al. (2009), that 
considered a slightly smaller black hole mass. Fig. [T| shows also 
the torus, disk and corona emission (unchanged for both states), 
the synchrotron component and the SSC and EC contributions for 
the Jan 13 SED. Particularly important is the UVOT spectrum in 
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Table 3. Parameters of the models shown in Fig.[T]and Fig.|2] a: set of parameters used in Ghisellini et al. (2010), where a slightly smaller M was adopted. 
Col. [1]: date; Col. [2], [3]: dissipation radius in units of 10^^ cm and in units of Rs; Col. [4]: black hole mass in solar masses; Col. [5]: size of the BLR in 
units of 10^^ cm; Col. [6]: power injected in the blob calculated in the comoving frame, in units of 10^^ erg s^^; Col. [7], [8]: accretion disk luminosity in 
units of 10*^ erg s^^ and in units of LEdd; Col. [9]: magnetic field in Gauss; Col. [10]: bulk Lorentz factor at iJdiss; Col. [11]: viewing angle in degrees; 
Col. [12], [13]: break and maximum random Lorentz factors of the injected electrons; Col. [14]: peak synchrotron frequency (Hz, rest frame); Col. [15] 
and [16]; slopes of the injected electron distribution [Q(7)]. 
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Figure 3. The energy density measured in the comoving frame of the blob 
as a function of the distance form the black hole, Rdms ■ We show separately 
the contribution of the magnetic field (Ub), of the radiation coming directly 
from the disk (U^^), and its X-ray corona (U'^),of the radiation produced by 
the BLR (C/gLp^) and by the torus (f/j'p^). The two vertical grey lines mark 
the values i?diss lor the two states of the source. The Jan 13, 2013 flare 
corresponds to dissipation occurring just beyond the BLR. For illustration, 
we also show U^j^^ for a "ring-like" BLR lying in the disk, shaped as a 
torus of radius i?BLR and cross sectional radius /?blr/10. 



Date 


logPr 


logPe 


logPc 


logPp 


23/12/2008" 


44.72 


44.56 


43.97 


45.92 


23/12/2008 


44.39 


44.55 


43.61 


45.62 


13/01/2013 


45.53 


44.54 


44.67 


45.25 



Table 4. Jet power in the forni of radiation, Poynting flux, bulk motion of 
emitting (therefore relativistic) electrons and cold protons, assuming one 
proton per emitting electron, a: set of powers derived in Ghisellini et al. 
(2010), where a slightly smaller black hole mass was adopted. 

Mgll lines. Given the uncertainties, these values are well consistent 
with ours. 

To reproduce the two states of the source, we have chosen to 
modify a minimum numbers of parameters. Besides Al and M, the 
models share the same magnetic field profile (i.e. the same Poynting 
flux Lfl = R^r^ see Tab.|4ll. The main changes concern the 
typical energies of the injected electrons, higher in the high states, 
the value of Rdias and the injected power. We have furthermore as- 
sumed that in the high state the buUi Lorentz factor is somewhat 
larger (r=16 vs. 13). As shown in Fig.|3] all these changes are con- 
sistent with the idea that the two states correspond to a different 
i?diss: for the "low" state, i?diss < -Rblr, the radiative cooling is 
strong, electrons cannot reach high energies and the overall SED is 
"red". By moving from 600 to 1600 Schwarzschild radii, the radia- 
tion energy density drops, letting the electrons to reach higher ener- 
gies. The magnetic field also decreases [B cx (-Rdissr)^^], but the 
change of the radiation energy density is more drastic, malsing the 
Compton dominance (ratio of the inverse Compton to synchrotron 
luminosities) to decrease in the high state. 



Dec 2008, that can be interpreted as the pealc of the accretion dislc 
component. The Shaw et al. (2009) optical spectrum is visible in 
Fig. (2] Its flux is larger than in Dec. 2009, there is a contamination 
from the synchrotron component at low frequencies, an upturn at 
large frequencies and the broad Mgll and H/3 lines are well visible. 
Their luminosities can be used to reconstruct the luminosity of the 
entire BLR (using the template of Francis et al. 1991 or Vander 
Berk et al. 2001) and then Ld- It agrees with the ones we have 
adopted, i.e. Ld ~ 3.5 x 10*^ erg s~\ Since the UVOT points 
show a peak in uF^, we find M — 7 x 10^ Mq. To show how well 
the black hole mass can be determined, we have plotted in Fig. |2] 
other two disk spectra, of the same Ld but with different masses 
(i.e. 0.4 and 1.4 billion solar masses). With these values the data 
are not well reproduced. By applying the virial method, Shaw et al. 
(2009) estimated M ~ (3 - 6) x 10* Mq, using the H/3 and the 



5 DISCUSSION 

Fig. m shows 7pcak as a function of the total radiation energy den- 
sity (magnetic plus radiative) as seen in the comoving frame of the 
source. The random Lorentz factor 7pcak corresponds to the en- 
ergy of the electrons emitting most of power (both in synchrotron 
and in inverse Compton). The grey points are the values of blazars 
analyzed by us in the past (see Ghisellini et al. 2012 and refer- 
ences therein), while the red squares corresponds to the 4 "blue" 
blazars discovered by Padovani, Gioinmi & Rau (2012). As these 
authors suggested, they are very likely FSRQs with a a very pow- 
erful optical synchrotron emission, that hides the emission lines. 
They should therefore coiTespond to the high state of PMN J2345- 
1555. The monitoring of these 4 blazars is too sparse to know if 
they are permanently "in the high state", or if we caught them by 
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Figure 4. Random Lorentz factor (7pcak) of 'he electrons emitting at tlie 
SED peaks as a function of the total energy density U' = U'^ + U'g as 
measured in the comoving frame. Small 7peak ™d large U' correspond to 
"red" blazars, while large 7pcak ^nd small U' correspond to "blue" blazars. 
The grey points correspond to blazar analyzed by us in the past (Ghisellini 
et al., 2012 and references therein), while red symbols refer to the 4 "blue" 
quasars studied in Padovani, Giommi & Rau (2012) and Ghisellini et al. 
(2012). PMN J2345-1555 is the first clear example of a blazar shifting from 
red to blue. 

chance. In any case, PMN J2345-1555 demonstrates that a single 
source can indeed vary not only its overall flux by orders of mag- 
nitudes, but also its overall "look" and "color". In our scheme this 
this is simply due to the change of the location of the dissipation re- 
gion by a factor less than 3 (from 600 to 1600 Schwarzschild radii 
in this case). This is sufficient to shift a blazar, in Fig.|4] from the 
red (7peak^lO^ and [/'>0.3 ) to the blue zone (7pcak?^10'^ and 
?7'<0.3 erg cm~'^). The top left region of the plane is populated 
mainly by low power, high energy peaked BL Lacs, lacking strong 
broad emission lines. For them, we do not expect strong variations 
of the cooling process when Rdiss changes (since there is no BLR). 
Variations are however possible due to changes in the acceleration 
process, in turn linked to the total power. Larger injected powers 
(and thus larger observed luminosities) could be associated to larger 
7pcak, as observed in Mkn 501 (Plan et al. 1998; Tavecchio et al. 
2001). 

The fact that PMN J2345-1555 is the first FSRQ "caught in 
the act" suggests that it is a rare phenomenon. Furthermore, out of 
a sample of hundreds of FSRQs detected by Fermi, only 4 were 
found to be blue by Padovani et al. (2012). This leads us to con- 
clude that only a few per cent of FSRQs dissipate most of their en- 
ergy beyond the BLR, either occasionally or permanently. In turn, 
this suggests that the region where the jet dissipates the most is in- 
deed well confined at i?disaSsl0^7?s- Larger radii (even by a small 
amount), corresponding to a blue FSRQ, are rare, although not im- 
possible. 

We envisage two possibilities to have a dissipation always well 
localized. The first is internal shocks: two shells initially separated 
by Ai? and moving with different F will collide at i?diss ~ AT?. 
If Ai? = a few Rs, and F ~ 15, then i?diss ~ IQ^Rs- Larger val- 



ues are possible for larger F or for larger AR. Alternatively, -Rdiss 
may be set by the pressure balance between the jet and the external 
medium. In this case 7?diaa is larger for a larger jet pressure, and 
then for a larger jet power 
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